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Abstract 



This note is motivated by the article by F. Lutscher, E. Pachepsky, and M. Lewis (2005), The Effect 
of Dispersal Patterns on Stream Populations SIAM Rev. Vol. 47 No. 4 pp. 749-772 on the drift paradox. 
We consider the case of a regime switching probabilistic model for a population of organisms living 
in a one dimensional environment with drift towards an absorbing boundary of the type introduced by 
C^ Lutscher et al. (2005). In particular, the two regimes consist of birth/death style demographics governing 

^^ the evolution of the immobile regime, and some form of advection-dispersion governing the evolution 

^> of the mobile regime, together with a regime-switching mechanism linking the two. In the present note 

n , it is shown for the regime-switching model, and in contrast to the results in the afore cited work, for any 

• finite advection speed, no matter how large, there is a finite critical length of the domain above which 

"^ the population can persist. 

1 , .,^«.. -«,.,..... 



J> In the paper Lutscher et al. ( 2005[ l, a benthos-drift model is presented, wherein a population of organisms 



OO existing in an environment subject to uni-directional flow is modeled as occupying two regimes: the benthic 

Z^ organisms are reproductive and non-mobile (these organisms are said to occupy regime in our terminol- 

7—i ogy), and the mobile organisms are non-reproductive (regime 1, here). Switching from the non-mobile to 

■^ the mobile regimes is taken to occur at some constant rate, here called Aq, and switching from the mobile to 

O the non-mobile regime also occurs with with constant rate, here called Ai ; such models fall within a general 

_■ class of so-called regime or hybrid switching models, see [Yin and Zhu|([2010|). 



It is argued in (Lutscher et al. ( 2005| ), p. 755) that, "if the rate at which individuals move times the prob- 



. ^H ability that they leave the domain during dispersal exceeds the population growth rate, then the population 

^\ will go extinct". In particular, this provides a critical equation for length as a function of the growth rate and 

j^ drift speed in the form 

r = Aovr(Ai, L,y) (1.1) 

where L ("length") is the location of an absorbing boundary, Y = {Yt : t > 0} is a stochastic process 
(to be described further) that describes the motion of individuals in the mobile regime, and the "wash-out 
probability"is given by 

7:{\i,L,Y) = P{Tl<Si), (1.2) 

where T^ is the first passage time to the absorbing boundary at L, and ^i, the regime switching time, is an 
exponentially distributed random variable with parameter Ai representing the amount of time an individual 



will remain motile before settling back into the benthos. In Lutscher et al. (2005 1 this formula (1.2i is 
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obtained by means of a spectral analysis of another model (1.3 1 for population density u that in turn is 
argued to yield ( |1.1[ ), where 

du(x,t) [^ 

— -^^ =ru{x,t) - Xu{x,t)+ Xk{y,x)u{y,t)dy, (1.3) 

for r, A > 0, and a transition probability kernel k{y, x). We consider the underlying regime switching model 
for the mobile and immobile populations from ( [Lutscher et al. ( |2005| ), p. 755) together with this formula 



( |1.1| ), in place of ( 1.3 1, as the starting point for an alternative analysis of the drift paradox. In particular, we 



will see that the critical length so-determined differs substantially from that determined by ( 1.3 1. 



The drift speed is determined by the following general assumption regarding the nature of the advection- 
dispersion in the mobile regime. 

The Drift Hypothesis (DH): Assume Y = {Yt : t > 0} is a stochastic process starting at Yq = and 
having independent, stationary increments with K{Yt-^.s — Ys) = vt < tJ < oo. Moreover, assume that 
the paths ofY are right-continuous without upward jumps. 

The stochastic process Y represents the motion of individuals in the mobile phase. Such processes are 
said to be spectrally negative Levy processes. The point of the spectral negativity assumption is to exclude 
motile particles that can jump over L without being absorbed there. We note that such a process is non- 
explosive under the assumption of right-continuity in (DH). 
The Levy-Ito decomposition for this process is 

Yt = vt + VDBt + Jt, (1.4) 

where i) > is deterministic drift to the right, D > 0, Bt is a standard Brownian motion, and Jt is a Levy 

process, which jumps only to the left. 

With V defined by (DH), it can be shown that 

v = v+ xn{dx), (1.5) 

Jo 

where n(-) is the Levy measure associated with this process. 

This quantity v is called the effective velocity, and it accounts for both deterministic drift to the right, and 

average Levy drift. Our assumption in (DH) that v > 0, amounts to a constraint on the pair v,n{-). 

The following simple lemma provides an alternative approach to the calculation of wash-out probabili- 
ties for this model in terms of the Laplace transform of the first passage time to the boundary evaluated at 
Al. 

Lemma 1.1. The wash-out probability is given by, 

7r(Ai,L,y) =Ee"^i^^. (1.6) 

Proof. 

7r(Ai, L,Y) = P{Tl < Si) = E{P{Tl < Si\Tl)) = E(e~^i^^). (1.7) 

D 



This is especially relevant when viewed in terms of the results obtained by Pakes ( |1996 1 that can be 



used to provide an explicit role of the parameters in connection with the wash-out probability in Lemma [TT] 
Namely, one has 



E {e-'^^) = e^'^'-'\ (1.8) 

where 6 = r]{s) is the unique solution to ip{6) = —s, and ip{9) is the cumulant function of the process Y: 

^(e)=iv+ xn{dx)]e + {1 - e-^"" - ex)n{dx) + (1 - e-^^)n(dx). (1.9) 
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In particular, using ( 1.8 1 it follows that 
Corollary 1.1. 



ETl = -, (1.10) 

V 



where v is the effective velocity, v = v + L xn{dx). 

Lemma 1.2. ThefunctionL — t- tt{Xi,L,Y), for a fixed value of Xi, is continuous and non-increasing, with 
7r(Ai, 0, Y) = 1 and limL^oo'^{Xi, L, Y) = 0. 

Proof To see that 7r(Ai,0,y) = 1, note that Yq = implies Tq = so that -k{Xi,0,v) = P(0 < 
^i) = 1. That TT is a non-increasing function of L follows easily from the fact that Li < L2 implies 
Tl, < Tl2 so that Tl^ < Si implies Tl, < Si and we have P{Tl2 < Si) < P{Tl, < 5i). To prove 
that /im/,_j.oo7r(Ai, L, Y) = 0, first note that T^ — )• 00 as L — )• 00 by the non-explosivity of the underlying 
Levy process Y under hypothesis (DH). So the asserted limit follows from Lemma [TTT] and the dominated 



convergence theorem. The continuity of vr as a function of L is also evident from the formulae of Pakes 



(19961. D 



In the case that < Ao < r, as noted in Lutscher et al. ( 2005 1, the population growth is obviously 



sufficient to prevent extinction without further need for analysis. Thus we restict attention to Aq > r and 
obtain the following main consequence of the previous lemmas. 

Theorem 1.1. Under (DH) with Aq , Ai , r > and r < Aq, there is an L = Lc such that r = Aovr(Ai , Lc, Y). 



Proof. The asssertion follows immediately from Lemma 1 .2 and the intermediate value theorem. D 



Example 1: Suppose that the process governing the mobile phase is Brownian motion with drift v and 
dispersion coefficient D > 0. Then we have 



,vL _ Vv'^ + 2DXi 
D W 



P{Tl < Si) = /tJAi) = exp(— - '^-^ ^L) (1.11) 



so that the critical equation ( |1.1| ) takes the form 



HyJ = {^ ^ )L, (1.12) 



so that 

ln(^)^ / 2DX 



\i 
In particular it follows that in the limit as f — )• 00 one obtains 



L = ^(yi + -^ + l).. (1.13) 



L^^ln{—)v. (1.14) 

Al r 



Observe that if D = then the asymptotic formula ( 1.14| ) is in fact the exact equation. So this example 
also covers the case of deterministic motion in the mobile phase but, more to the point, also shows that large 
enough v nuliifies the dispersion effect. 

We see here that, although these results are based on both the underlying regime switching model intro- 



duced in (Lutscher et al. (2005 1, Sec 3.1), and the characterization of the critical length given by expression 
(3.8) there, the conclusions above differ qualitatively from Lutscher et al. ( |2005j ) concerning equation ( |1.1| ) 
when calculated in the context of ( |1.3| ). Specifically, for any finite advection velocity, no matter how large, 
there will be a length that satisfies ( |1.1[ ), whereas in Lutscher et al. (2005) it was found that there was a 
critical advection velocity Vc such that L ^ oo as v ^ v~, and above which no length will satisfy ( |l.l[ ); 
see Figure[T]for comparison to (Lutscher et al. (2005 1, Figure 1). 




0.5 1.0 1.5 

Effective Velocity, u 



Figure 1 : Critical length as a function of drift velocity 



2 Some Concluding Remarks 

While the main result of this paper is somewhat surprisingly at odds with the formula (4.6) for the critical 



length obtained by Lutscher et al. ( 2005 1 as a function of drift speed, the resolution appears to be that the 



the regime-switching model and the model ( 1.3 1 are in fact quite different. 

In a still different but somewhat related framework to that oflLutscher et al.l (120051), IKestenl (11978) ob- 



tained the existence of a finite critical drift speed for a system of branching Brownian motions with negative 
drift and absorbed at the origin. While such results capture the spirit of the result in |Lutscher et al. ( |2005 1, 
this involves non-linear production of offspring but without regime- switching. The following simple result 



is provided to simply clarify the respective roles of linear and non-linear production from a probabilistic 
perspective. The case of non-linear production goes back to Henry McKean's probabilistic representation 
of solutions to the Kolmogorov-Petrovski-Piskunov-Fisher (KPPF) reaction-diffusion equatioi|^in |McKean] 
( [1975[ ). While we expect that the description given below for the linear case is known, we are unable to find 
a referenceo 

Theorem 2.1. Let Arf{x) = 2Df"{x) + vf'{x) + rf{x) be the infintesimal generator of the transition 
operator T[, where r > 0, D > 0, ij G M. Also let Y be a Markov process with infinitesimal generator 
A{y = ^D-^ + v-g^ and transition operator T^, and let N = {Nt : t > 0} be a Yule-Furry branching 
process with rate r > 0, which is independent ofY. Then for f G D^ one has that Mt = Ntf(Y{t)) — 
Jq ArNsf{Ys)ds, t > 0, is a martingale. 

Proof. Let Tt = (t{Ns, Ys : s < t). Then we have, for all < s < t, 

E(Mt|J-,) = ¥.{Nt\Ts)nf{yt)\^s)- r Ns>fiYs')ds'- [ EiNs:Arf{Ys>)\Ts)ds' 

Jo Js 

rs ft QJ^O f(Yi) 

= N,e^^'-'^Tlj{Y,)- j ArN,,f{Ys>)ds'- J e^(^'-^)iV,(^^^^ 

+rT^,_J{Y,))ds' 

= Nse'^'-'^TlJiYs) - r ArN,f{Ys>)ds' - {e<'-'^'^NsTlj{Ys) - Nsf{Y,)) 

Jo 

= NJ{Ys) - r A,N,,f{Y,,)ds' = Ms. (2.1) 

D 



It follows from Theorem [zT] that letting y(j) =Y, j = I, . . . ,Nt, 

Nt 

u{x,t) = E J^/(y/^')) = E,.iVt/(yt), t > 0, (2.2) 

j=l 

solves 

— = Aru{x, t), u{x, 0+) = f{x), X e R, t > 0. (2.3) 

In the case of KPPF the branching offspring evolve independently along paths distributed as their parent, 
whereas in the case of linear production the offspring particles are clones following precisely the same path 
as their parent. This is an essential difference between the stochastic branching and (linearized) dynamic 
growth models. In the case of branching Brownian motions with negative drift and in the absence of regime 



switching, Kesten ( 1978 1 shows that at sufficiently high drift rate all particles will eventually be absorbed 



at the outlet regardless of L, namely \v\ > \/2Dr, while there is a chance to persist at lower drift speeds 
toward the absorbing boundary. For the linear production rates analyzed here and in Lutscher et al.| ( [2005] ). 



because of regime switching and cloning the eventual absorption is sure to occur for any positive drift speed 
only if L is too small, i.e., it is always possible to find a (finite) length L large enough for persistence under 
arbitrarily specified drift and regime switching rates. 



'The KPPF equation is a reaction-diffusion independently introduced by Kolmogorov, Petrovski and Piskunov on one hand, 
and by Fisher on the other. It is sometimes simply referred to as KPP in the older mathematics literature. However Fisher's 
considerations in populaton genetics lead to the same equation, and more modern references either refer to FKPP or KPPF. 

^Of course, in the case of a negative rate in the linear case, the probabilistic representation is also well known in terms of 
'particle killing'via the Feynman-Kac formula. 
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